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ABSTRACT 

In heavy quark production at large Feynman x there are two hardness scales, 

Qh! one given by the heavy quark pair mass A^^ and the other by Aq^£,/(1 — x). When 

these two scales are comparable, the twist expansion of Perturbative QCD breaks 

K^ '. down. We discuss the dynamics in this new QCD limit, where yU^ = 7V1^(1 — x) 

. 5t ! is held fixed as 7W^ -^ oo. New diagrams are found to contribute, which can 

enhance the cross section above that expected for leading twist. The heavy quarks 

are produced by a peripheral scattering on the target of hardness /j?. This leads, 

in particular, to a nuclear target dependence of A^'^ at small /i. Qualitatively, 

the dynamics in the new limit agrees with earlier phenomenological models of 

"intrinsic" heavy quark production. 



• Lectures presented at the XXXII Cracow School of Theoretical Physics, Zakopane 1992. 



1. Introduction 

The application of Perturbative QCD (PQCD) to hard processes involving 
hadrons and leptons has been very successful. The predictive power of PQCD 
is based on the factorization theorem , according to which an observable cross- 
section a can be factorized into a product of universal structure and fragmentation 
functions times the hard constituent cross-section a. Generically, 

a = FaFba{ab^c + X)Dc (1.1) 

where Fa, Fi, are the single parton structure functions of the incoming hadrons 
("probabilities for finding the quarks/gluons a, b in the projectile and target") and 
Dc describes the hadronization of the produced quark(s)/gluon(s) c into the ob- 
served final state {e.g., jets or hadrons at large pj_). The structure and fragmenta- 
tion functions F, D are not calculable in PQCD, but are predicted to be universal, 
i.e., independent of the specific hard collision a. The "higher twist" corrections 
to the "leading twist" PQCD prediction (1) are suppressed by powers of the large 
momentum scale Q^, which characterizes the hard subprocess ab -^ c + X. 

Here I would like to discuss modifications to Eq. (1.1) that are expected, and 
observed, in the case that one of the partons a,b ot c carries a large fraction x of 
the available longitudinal momentum. In the x — *> 1 limit there is a new hard scale 
Aq(-.£)/(1 — x), and the corrections to eq. (1.1) are of order Aq(jjj/Q'^(1 — x). In 
the combined limit 

> with /i^ = Q^{1 - x) fixed (1.2) 

x^l ) 

2 

the twist expansion in fact breaks down . The higher twist contributions are of 
(9(l//i^), and hence not suppressed by an asymptotically large variable in this 



new QCD limit. Since large momentum transfers are involved, PQCD can still be 
used to analyze the process. As I shall discuss below, in this limit the cross-section 
cannot be expressed in terms of single-parton structure functions; it involves multi- 
parton distributions. This is similar to the case of exclusive {x = 1) hard scattering, 
which depends on the longitudinal momentum distributions of all (valence) quarks, 

3 
constrained to be at small transverse distances from each other . 

2. Why the a; ^ 1 Limit is Hard 

Consider the wave function 0(yp, n_|_) of the hadron h in Fig. 1, in a frame 
where h has a large longitudinal momentum p. We shall assume that describes 
the soft, non-perturbative part of the quark distributions, and is suppressed in 
the limits y -^ 0,1 and also for n_|_ -^ oo. The perturbative tail of the full wave 
function when the fractional momentum of one constituent x — > 0, 1 or when its 
transverse momentum j9_|_ — » oo is then generated by gluon exchange as shown in 
Fig. 1. 
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It is in fact straightforward to see that Fock states where one constituent 
carries a large momentum fraction a; ~ 1 must have a short life-time, and hence 
are calculable in PQCD. The (kinetic) energy difference between the final Fock 
state of Fig. 1 and the hadron is 
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where we assumed p » p±/{l — x). By the uncertainty principle the "life-time" 



* The following simplified discussion ignores multiple gluon exchanges, which will bring log- 
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arithmic corrections, in analogy to the treatment in Ref. . 



Tqq of the Fock state is then proportional to 1 — x: 



Tg,---^- ^^^ (l-x) (2.2) 



Furthermore, since the hfe-time of the x ^ 1 Fock state is brief, we would 
expect that the transverse distance r±_ in Fig. 1, between the quarks before/during 
the gluon exchange must be similarly short, so that the duration of the exchange is 
no longer than the life-time Tqq. This is readily verified because the gluon exchange 
amplitude depends on the constituent transverse momentum n±_ only through the 
energy difference associated with the dashed line in Fig. 1, 



2pAE ^ !^i±^ - !^i±^ - i^^^lll (2.3) 

1 — y 1 — X X — y 



For a; ^ 1, AE is independent of n_|_ for 

nl<0{pl/{l-x)) (2.4) 

Since the non-perturbative wave function of the hadron will cut off the integration 
over n_|_ before the limit (2.4) is reached, the nj_-integration can be factorized, 

dn±(j){yp, n_L) = 0(i/p, r_L = 0) (2.5) 

showing that the only hadronic Fock states which can generate the x -^ 1 per- 
turbative tail are those with a short transverse distance r^ oc 1 — a; between the 
quarks. Conversely, according to Eq. (2.4) the x — > 1 Fock states involve large 
transverse momenta n^ oc 1/(1 — x). 



It is important to note that the x ^ 1 Fock state, which is compact at its 
moment of creation, nevertheless quickly expands in the transverse direction. Dur- 
ing the effective life-time (2.2) of the Fock state, the "slow" quark can move a 
transverse distance 



P± 2px _^) r^ 2p_La^ 

p{l — x) m^ + p]_ m'^+p']_ 



R±-V^Tgq = — - ^ ^ ^ (I-X)- 0,0 (2.6) 



which for p_|_ = O^Aqcd) is of the order of 1 fm. This observation will be important 
in the following. 

Because of the hard scale (2.4), PQCD can be used to calculate the behavior 
of hadronic structure functions in the x ^ 1 limit. The result is given by the 
"spectator counting rules" 

^ oc (1 - x)2"-i+2|AA| (2.7) 

ax 

where Ug is the number of "spectator" partons (whose Xs — * 0) and AA is the 
helicity flip between the initial hadron and the observed quark (or gluon) carrying 
the large momentum fraction x. 

3. Breakdown of the Twist Expansion for x ^^ 1 

Consider the well-known process of Deep Inelastic Scattering (DIS) of leptons 
on hadrons. The diagrams are classified as "Leading Twist" (Fig. 2a) or "Higher 
Twist" (HT, Fig. 2b) according to whether the spectator partons in the target are, 
or are not, connected to the active (hit) quark (or to partons radiated from this 
active quark). 



In the Q ^^ oo (fixed x) limit, the HT diagrams are suppressed by 1/Q ■ 
This can be intuitively understood as follows. The hard Iq -^ d'q subprocess has 
a duration r ~ 1/Q- Any interaction with spectators as in Fig. 2b can affect the 
hard scattering probability only if it occurs within this short time interval r (later 
interactions will only modify the momentum distribution of the struck quark in 
the final state). But an interaction within a time-scale r is possible only if there 
are spectators within a transverse distance r^ r^ t r^ ^/Q- The probability for this 
is proportional to the transverse area vrrj^ ~ ^IQ\ which explains the suppression 
of the higher twist contributions. 

The above argument for the suppression of HT terms breaks down in the high 
X limit. As we argued earlier, the x ^ 1 Fock states in Fig. 1 are produced from 
compact hadron configurations, with a typical distance r^ ~ (1 — x)/ ji^ between 
the valence quarks, where /x ~ 1 fm^ . Hence in DIS with Q^ ~ f^'^/i^ ~^)i '^•C-5 i^^ 
the limit (1.2), the scale of the hard photon interaction is commensurate with the 
size of the valence Fock state, and the scattering is coherent over several quarks. 
This means that the DIS cross-section cannot be expressed in terms of single-parton 
structure functions in the limit (1.2), and the usual factorization (1.1) fails. 

The magnitude of the higher twist terms has been studied experimentally in 

7 
DIS as a function of a;. The results show that the HT corrections are important for 

X ^ 0.5, and have an x-dependence which is consistent with (1 —x)^^, as suggested 

15 8 

by the above qualitative arguments and explicit calculations ' ' . Similarly, in high 

9 10 
mass lepton pair production there is experimental evidence ' for corrections to 

the Drell-Yan process, which are in qualitative agreement with the expectations 

for higher twist effects . 



4. Dynamics in the new QCD Limit 
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We have studied the production of quark pairs with large invariant mass A1 
in the high x hmit corresponding to (1.2), i.e., for 

with /i^ = M^il - x) fixed (4.i; 

X ^ 1 



where x is the momentum fraction of the qq system. An exphcit calculation of 
all relevant diagrams was done for scalar QED. Here I would like to discuss the 
qualitative conclusions, which apply equally to QCD. 

Consider first heavy quark production in the standard QCD limit, 

M"^ ^ oo at fixed x. (4.2) 

The usual lowest order diagram describing the fusion process GG -^ cc is shown 
in Fig. 3a. The virtualities of both gluons range up to O^Ai"^) - their momentum 
distributions are given by the projectile and target gluon structure functions eval- 
uated at a scale Q^ = A4^. Similarly in Fig. 3b, the higher order (in as) process 
qG -^ qQQ is given by the antiquark and gluon structure functions of the projec- 
tile and target, respectively, evaluated at a scale A4^. Both diagrams in Fig. 3 
are leading twist - they involve only one of the partons (a gluon and an antiquark, 
respectively) in the projectile, and one gluon in the target. 

The dynamics of the new limit (4.1) differs in several respects from the above. 
There are two types of leading order diagrams, the "extrinsic" and "intrinsic" ones 
shown in Figs. 4a and 4b, respectively. In extrinsic diagrams the heavy quark 
pair couples to only one parton in the projectile, while in intrinsic diagrams it 



couples to several. Since the produced QQ pair carries almost all of the momentum 
in the final state {x -^ 1), the light valence quarks g, g are effectively stopped. 
The light quarks then give a big contribution to the energy of the intermediate 
states, 2pAE ~ Pj_/(1 — x) (cf. Eq. (2.3)). The production cross section is 
dominated by values of the light quark transverse momentum pj_ where this light 
quark contribution to the energy difference is of the same order as that of the heavy 
pair, 2pAE ~ M^. Hence (for fi > J^qcd) 

pi^ M'^{l-x)=fi^ (4.3) 

Because of their small fraction ~ {1 — x) of the projectile momentum, the stopped 
light quarks can move a considerable distance -R_l in the transverse direction, even 
during the brief life-time of the virtual qqQQ Fock state. According to Eq. (2.6), 
i?_i_ ~ l/p_L ~ 1//^- Now a Fock state of this transverse size can be "resolved" 
by a target gluon of transverse momentum £_|_ (or virtuality £^ ~ ~^j_) of order 
£_i_ ~ l//i. Hence the hardness of the scattering from the target does not increase 
with the heavy quark pair mass M.. To be able to describe the target scattering 
using PQCD, we would have to choose fi >> Aqcd in the limit (4.1). In general, 
however, we must conclude that arbitrarily heavy quarks can be, and are, produced 
at high X by soft peripheral scattering. Such soft scattering is surface-dominated 
for nuclear targets, i.e., a oc A"^'^ as observed for J/^ production at large x (see 
below). 

Note that soft scattering is kinematically allowed only at sufficiently high en- 
ergies. The minimum longitudinal momentum transfer from a stationary target 



* Equivalently, we could say that an interaction of this hardness with the target deflects the 
stopped light quarks sufficiently to break up the Fock state and materialize the heavy quark 
pair. 



required to put a heavy quark pair of mass A4 on its mass shell is of order A4 /2p, 
where p is the laboratory momentum of the projectile. For charm production the 
minimum momentum transfer is below 50 MeV already for p > 100 GeV. 

The large transverse size R± of the light quark distribution also explains why 
the scattering dominantly occurs off the light quarks, as shown in Fig. 4. The 
heavy QQ pair has a small transverse size hj_ ~ 1/7W. A target gluon can couple 
to, and resolve , the QQ pair only provided it has a commensurate wavelength, 
i.e., i± ^ Ai as indicated in Fig. 3a. This is much larger than the i± ^ fi required 
to resolve the light quarks. Hence the gluon fusion diagram of Fig. 3a, which 
gives the leading contribution in the usual, fixed x, QCD limit (4.2), is actually 
suppressed by l/7Vf compared to the diagrams of Fig. 4 in the new limit (4.1). 

The Fock state of the projectile hadron from which the heavy pair is produced 
must have a small transverse size r^ ~ (1 — x)/ jj? ~ 1/ M.'^ (cf. Fig. 4). The 
argument is the same as the one already given in Section 2, leading to the estimate 
(2.4) of the valence quark transverse momentum before the virtual creation of the 
heavy quark pair. 

The extrinsic (Fig. 4a) and intrinsic (Fig. 4b) diagrams are of the same order 
in as and have the same behavior in the limit (4.1). However, we note some 
qualitative distinctions between these two classes of diagrams: 

(i) The intrinsic diagrams do not contribute significantly to lepton pair pro- 
duction, since two photon exchanges would be required. This may explain the very 
different A-dependence of lepton pairs, as compared to the J/ip (see below). 

(ii) The intrinsic diagrams give rise to a non-trivial phase in the leading order 



* If the virtual QQ pair is in a color octet configuration it will also interact coherently with 
soft gluons. In this case the virtual pair behaves like pointlike gluon, and is not materialized. 



amplitude. As shown in Fig. 5, in an intrinsic diagram the QQ production can 
proceed through two consecutive real processes. First an on-shell QQ pair with 
X < 1 is formed, and later the pair is accelerated to a; ^ 1 via an interaction 
with the second valence quark. This dynamical phase can be of importance in 

12 13 

polarization phenomena, which have been observed ' to be enhanced at large x. 
(iii) The extrinsic diagrams dominate over the intrinsic ones both in the exclu- 
sive (/i <^ ^QCD in Eq- (4.1)) and inclusive (/i ^ Aqcd) limits. In the extrinsic 
diagram of Fig. 4a, the virtuality of the gluon exchange between the light valence 
quarks is set by 1/t^ ~ Aq(jjj/(1 — x), while that of the gluon connecting to the 
heavy quark pair is set by 7W^. In the intrinsic diagram of Fig. 4b, on the other 
hand, both gluons have virtualities determined by the larger of these scales. In the 
scalar QED model calculation of Ref. 2, the intrinsic diagrams nevertheless domi- 
nated at intermediate values of /i^. It would clearly be important to determine the 
magnitude of the intrinsic contribution in QCD. 

5. The Production of Charmed Hadrons at Large x 

I _ 14 

The data on charm production in e^e annihilations and in photoproduc- 
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tion agrees well with PQCD at leading twist. The charm fragmentation function 
D{c -^ D + X) can thus be determined from these reactions using the factorization 

16 
formula (1.1). The e+e^ data show that the charmed hadron carries an average 

fraction (z) ~ 70% of the charm quark momentum. The momentum distribution 

17 

is often parametrized in terms of the "Peterson" function , 

""^/^^'^ " Z{1 - 1/z -\e/(l - Z)r ^'-'^ 

where ec — 0.06. According to the QCD factorization theorem, this fragmentation 

10 



function should be the same in all hard processes, regardless of how the charm 
quark is produced. 

There are several experiments on hadroproduction of charm ' ' ' ' which 
give much larger cross sections at high x than expected from leading twist QCD, 
i.e., from Eq. (1.1) with a fragmentation function of the form (5.1). Recently, 

23 24 

data with good statistics on n^A ^ D + X ' clearly shows that Eq. (1.1) 
underestimates the charm cross section already at medium values of x ^ 0.2. As 
shown in Fig. 6a, the experimental D-meson distribution actually is similar to 
that predicted by Eq. (1.1) for the charm quark ' . Thus the fragmentation 
function D{c ^ D + X) must be assumed to be close to 5(1 — z) in order to get 
agreement between experiment and theory in charm hadroproduction at medium 
X. Hadroproduced charm quarks appear to fragment in a strikingly different way 
from what they do in e^e^ annihilations, where the fragmentation function is to a 
good approximation given by Eq. (5.1). This implies a breakdown of the leading 
twist factorization formula (1.1). 

The different form of the fragmentation function in e^e~ induced, compared to 
hadron induced, charm production has a natural explanation ' . In hadroproduc- 
tion, the charm quark can coalesce with light spectator quarks from the projectile, 
which move in the same direction and with similar velocity as the charm quark. 
The charmed hadron formed this way will have the same velocity as the charm 
quark, resulting in a fragmentation function close to the (5-function suggested by 
the data. 

The existence of coalescence between charm quarks and light valence quarks 
from the projectile is suggested also by the "leading particle effect" . Those charmed 
hadrons that have a valence quark in common with the projectile are experimen- 

11 
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tally found ' ' to have a harder x-distribution. In Fig. 6b, c we show the size of 

23 25 

this effect as measured by E769 , and as expected in a model with coalescence . 
From this model one can also see that the E769 data is not very sensitive to the 
estimated contribution of intrinsic charm, which is important only at the highest 
values of x. 



6. The A Dependence of Quarkonium Production 

The coalescence of heavy quarks with light comovers has an indirect effect on 
the production of quarkonium states, such as the J /ip and T. When many co- 
moving light quarks are present, as in the fragmentation region of heavy nuclei, 
the heavy quarks may preferentially coalesce with comovers rather than bind to 
each other . This leads to a suppression of quarkonium production in nuclear 
fragmentation regions, which has been observed for the J /ip in central heavy ion 
collisions . A similar suppression is also observed in "backward" (x < 0) pro- 
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duction of both the J/il) and the T in pA collisions , for which the alternative 
explanation in terms of a quark gluon plasma seems unlikely. 

The forward {x > 0) production of the J/ip has been measured with high 
statistics for both pion and proton beams on a variety of nuclear targets. ' ' In 
this case the effect of coalescence should be small, since the cc state is produced 

in the fragmentation region of a hadron, and has relatively few comovers. This 

34 
data nevertheless gives direct evidence for the breakdown of the leading twist 

approximation at large x. In the factorized formula (1.1), the nuclear target A- 

dependence can only appear through the target structure function F^/^(a;2), where 

X2 is the momentum fraction carried by the target parton b. Ratios of J/ip cross 

sections at different production energies but at the same X2 should therefore be 

12 
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independent of the nuclear number A. As shown in Fig. 7 for the NA3 data 
on n^A -^ J/ip + X at 150 and 280 GeV, the ratio of cross sections on Hydrogen 
and Platinum does not in fact scale as a function of X2- Thus the leading twist 
factorization (1.1) fails. A similar result was obtained by combining pA data from 
NA3 and E772^^ 

The failure of factorization occurs at the smallest values of X2, i-e., for large 
Feynman x of the J/iJ, since X2 — Mj,^/xs. In this region of x the nuclear target 
dependence of the J/ip cross section also is not linear in A, as expected for hard 
QCD processes which occur incoherently off all partons in the nucleus. If the J/ip 
cross section is parametrized as 

CThA = o-/iiV^" (6.1) 



one finds ' ' that a = 0.7 . . . 0.8. 

35 
At small values of X2, one does expect a ^ 1 due to parton shadowing , as 

36 
observed in deep inelastic lepton scattering (DIS) and in lepton pair production 

(DY). However, shadowing appears to be a Q^-independent, leading twist effect 
associated with the nuclear structure function. It thus cannot account for the 
breakdown of factorization observed in J/ip production. Moreover, the shadow- 
ing effect seen in DIS and in DY is a fairly small, 10 . . . 30% effect, whereas the 
suppression of J/ip production amounts to a factor ~ 3 for large nuclei. 

It was recently suggested ' ' that the nuclear suppression of J/ip production 
could be due to energy loss of the incoming and outgoing partons while propagat- 
ing through the nucleus. However, very little energy loss is observed for the struck 

40 41 

quark in the DIS process , as well as for the incoming quark in the DY reaction . 
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This shows that high energy quarks suffer only insignificant energy loss in the nu- 
cleus. Actually, this fact is a direct consequence of the uncertainty principle, which 
forbids enhanced energy loss from multiple collisions occurring within the forma- 

42 

tion zone of the radiated gluons . Repeated radiation from collisions separated 
by a time interval At is allowed only provided AEAt ^ 1, i.e., when the energy 
difference resulting from the gluon radiation AE is big enough for the multiple 
scatterings to be resolved. In a nucleus. At ^ Ra, where Ry\ is the nuclear radius. 
The emission of a gluon with transverse momentum p±_ and energy Eg results in 
AE ~ p^/2Eg. Hence normal, soft collisions in the nucleus with (p]_) ~ 0.1 GeV 
can only lead to a finite energy loss in the laboratory frame, 

Eg^^ (pI) Ra^I-5 GeV. (6.2) 

For high energy partons this fixed energy loss is not significant, and it cannot 
explain the nuclear suppression of J/ip production. 

When the fixed energy loss (6.2) is small compared to the energy of the J/ip, 
we expect to see Feynman scaling (in x) of the J/ip cross section. This is in- 
deed observed in the data (c/. Fig. 7). Hence it is most natural to discuss the 
J/ip x-distribution from the point of view of the projectile wave function. This 
brings us back to our earlier discussion of heavy quark production at large x. The 
74-dependence of the J/ip data is indeed one of the strongest arguments for the 
relevance of the new QCD limit (4.1), where A^^(l — x) is held fixed as A^^ -^ oo 
and X -^ 1. 

In Section 4 we saw that a virtual heavy quark pair with x near 1 can be put 
on its mass shell (produced) by a soft gluon scattering off the stopped light quarks. 
If 1 — X = jj? /Ai"^, the light quark distribution has a transverse size l//i, and the 

14 



hardness of the target interaction is Q = /i (c/. Eq. (4.3)). Hence for small 
/i ~ Aqcd, the scattering will be surface-dominated in a big nucleus, and the cross 
section gets an A ' nuclear target dependence. There should thus be a smooth 
transition from an A behavior at small x (hence large /i) to an A ' behavior at 

43 

large x, in agreement with the data. 

31 

The NA3 data was analyzed assuming the existence of two components in 
the cross section, a "hard" component with a = .972 in (6.1) and a "diffractive" 
component with a = .77 (for the 7t~A data) or « = .71 (for the pA data). A 
quantitative fit showed that the hard component was in good agreement with ex- 
pectations from the GG -^ cc and qq -^ cc fusion processes. Hence the "diffractive" 
component, which dominates for x ^ 0.6, appears to be an excess over the leading 
twist QCD fusion contribution. Note that explanations of the anomalous nuclear 
dependence in terms of energy loss or breakup of the J/ip would not lead to any 
excess in the large x production cross section on light targets. 

Because the target interaction in the x ^ 1 limit (4.1) can be soft (£_|_ ~ /i in 
Fig. 4), the average transverse momentum (A;_|_) of the J/ip should decrease from 
(A;_i_) = 0{mc) at low a; to (A;_|_) = 0{fi) at large x. This effect has been seen in the 

12 

J/ip data . Furthermore, the anomalous A-dependence, i.e., a < 1 in Eq. (6.1), 
31 32 

is observed ' only at low k±. A similar decrease of (A;_l) with x is expected 
also in open charm production. In addition, coalescence between heavy and light 
quarks occurs mainly for heavy quarks with low k±, comparable to the light quark 
transverse momenta. This also tends to decrease the (A;_|_) of charmed hadrons at 
large x. 

For T production, the available data is not at large enough x to make /i^ = 
(1 — x)A4'^ small. Hence the corrections to the leading twist results should not be 
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very significant. The A-dependence of tlie T cross section is in fact found to be 
much closer to A^ than is the case for the J/^. 



7. Concluding Remarks 

We have discussed the corrections to the twist expansion of PQCD for hard 
processes in the kinematic region where some of the constituents (or hadrons) carry 
a large fraction x ^ 1 of the available energy. We found that the limit (1.2), where 
the hard scale Q^ (or Ai"^) is comparable to the scale Aq^j~,/{1 —x), is particularly 
interesting. The production mechanism of a heavy quark pair of mass Ai becomes 
peripheral in this limit, with the hardness /i^ of the target interaction being given 
by yU^ = 7Vf^(l — x). For fi ^ Aqcd this means that the interaction in the target 
is soft. In this situation the factorization formula (1.1), according to which the 
scattering occurs incoherently off single partons in the target, breaks down. In 
particular, it implies scattering from the surface of large nuclei, and a consequent 
A"^'^ nuclear target dependence. 

On the other hand, only very compact Fock states of the projectile, where the 
transverse separation of the light valence quarks is of order l/TW, can participate 
in the production. For large x several valence partons of the projectile are involved 
in the heavy quark production process, again breaking the leading twist, single 
parton factorization (1.1). In this respect the dynamics is similar to that in the 
exclusive {x = 1) limit. The hard production vertex can be calculated in PQCD, 
given the longitudinal momentum distributions of the quarks in the transversally 
compact Fock states of the projectile. So far, this calculation has only been carried 

2 

out in a scalar QED model, however. 
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Phenomenologically, significant deviations from tlie leading twist approxima- 

7 

tion liave been seen at large x in deep inelastic lepton scattering , in lepton pair 

q 1 Q Ql QO 'V\ 

production ' , in J /tp hadroproduction ' ' and in open charm hadroproduc- 
tion ' ' ' ' . Recently, firm evidence was obtained ' that there are significant 
corrections to the leading twist PQCD cross section for open charm hadroproduc- 
tion even at moderate x ^ 0.2. The most natural explanation for this effect is the 
coalescence of the charm quark with co-moving light spectator quarks ' . Coa- 
lescence is a soft process, which is not directly calculable in PQCD. It allows the 
charm quark to maintain its velocity during hadronization, and can thus effectively 
be described by a fragmentation function D{z) ~ 5(1 — z). 



A significant change in the velocity of a heavy quark always requires large 
momentum transfer - hence this is a process that is associated with the hard pro- 
duction vertex, not with soft hadronization. Some fragmentation schemes based on 
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the string model ' give rise to charmed hadron distributions which are consider- 

22 24 

ably harder than the distribution of the charm quark ' . Such mechanisms seem 
to go beyond the realm of soft physics for which the models were intended, and are 
thus unreliable. The methods that we have described above illustrate how the cross 
section of heavy quark production at the largest values of x can be calculated in 
QCD. While the size of the new, intrinsic contributions to charm production has 

25 43 

not yet been determined from the theory, phenomenological estimates ' show 
that they are important for x ^ 0.5. Several experiments have in fact reported 
larger charm hadroproduction cross sections at high x than would be expected 
from the leading twist approximation (1.1). These data await confirmation by 
upcoming, high statistics experiments. 
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FIGURE CAPTIONS 

Fig.l The x — > 1 limit of a hadron structure function is generated by perturbative 
gluon exchange. The transverse size r±_ of the initial Fock state is small, r^ ~ 
(1 — x)/pj_. The virtual state with large x has a larger size R± ~ l/p±, due 
to the transverse motion of the stopped quark carrying the small momentum 
fraction 1 — x. 

Fig. 2 Deep inelastic lepton-hadron scattering, (a) A leading twist diagram, (b) A 
higher twist diagram. 

Fig. 3 Leading twist diagrams in heavy quark production at fixed momentum frac- 
tion X of the heavy pair, (a) A lowest order GG -^ cc diagram, (b) A higher 
order qG — > qQQ diagram. The gluons have virtualities ranging up to the 
mass scale Ai of the heavy pair. There is no restriction on the transverse 
size of the initial Fock state. 

Fig. 4 Leading order diagrams in the limit (4.1). In the extrinsic diagram (a) the 
produced heavy quark pair couples directly to only one parton in the projec- 
tile. Diagram (b) is intrinsic, as the QQ pair couples to two partons. In the 
standard limit (4.2), both diagrams (a) and (b) would be classified as higher 
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twist, and would be suppressed by l/TW . The gluon from the target has 
a virtuahty of 0{fi), while the two other gluons have virtualities of 0{Ai). 
The transverse sizes r_|_ and R± of the initial and final Fock states are as in 
Fig. 1, and the size of the QQ pair is h± ~ 1/A^. 

Fig. 5 In an intrinsic process the intermediate state indicated by the vertical dashed 
line can be on-shell, resulting in an imaginary part of the amplitude at leading 
order. 
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Fig. 6 xp distributions of D mesons produced in tc^A collisions at 250 GeV/c 
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(E769 data ). (a) The solid line shows the prediction of leading twist QCD, 

17 
using a Peterson fragmentation scheme for c ^ D + X which fits the data 

on D production in e^e^ annihilations. The dashed curve shows the effect of 

adding an estimated intrinsic charm contribution to the fusion cross section. 

The dot-dashed curve results from fusion + intrinsic charm at the quark 

level, i.e., for 5-function fragmentation, (b) D^ production, for which the 

curves include both the fusion and intrinsic production mechanisms. The 

solid curve is obtained with Peterson fragmentation and the dashed curve 

with (5- function fragmentation, (c) As in (b), for D^ production. 

Fig. 7 The ratio R = Aa{pp -^ J/ip + X)/o'{pA -^ J/ip + X) of inclusive J/ip 

31 
production cross sections on Hydrogen and Platinum . In (a) the ratio is 

plotted as a function of the Feynman xp of the J/ip, and in (b) as a function 

of the momentum fraction X2 of the target parton . 
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